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We review recent results obtained for a class of one-dimensional quantum models 
with 1/r 2 long-range interaction. Based on the asymptotic Bethe-ansatz solution and 
conformal field theory, we study critical properties of the continuum boson model, the 
SU(z/) spin chain, the OSp(V,l) supersymmetric t-J model, and a new hierarchy of 
models related to the fractional quantum Hall effect. We further investigate the class 
of 1/r 2 models with harmonic confinement by means of a newly proposed method of 
the renormalized-harmonic oscillator solution. 



PACS numbers: 



Typeset Using REVTEX 



1 



§1. Introduction 



One-dimensional (ID) many-body systems with long-range interaction of 1/r 2 type 
||r)H3li have been studied extensively in connection with fundamental notions in condensed 
matter and statistical physics. Even in classical models [j4])-|7l, one can find remarkable 
results. For example, a path-integral formulation of the Kondo problem results in the fer- 
romagnetic Ising model with 1/r 2 interaction Also, the antiferromagnetic Ising model 



shows a typical example of the devil staircase |6j 

Quantum models with 1/r 2 interaction have attracted renewed interest considerably [1-3, 
8-42], revealing new interesting aspects of integrable systems, and thus providing us with a 
paradigm of fundamental ideas actively discussed in condensed matter physics, such as the 
random matrix pJpJH , the Gutzwiller state for correlated electrons [|3j|- |10)|| , the fractional 
quantum Hall effect (FQHE) | |17)| - |22)| ], the level statistics for disordered systems [ 23j[ -f25)|, 
etc. 

A common feature in quantum 1/r 2 systems is that the ground state is exactly given 
by a Jastrow-type wavefunction [1-3, 8-22], namely a product of two-body functions. This 
characteristic nature of the wavefunction should be closely related to the integrability of 
1/r 2 models because it implies that the two-body scattering is essential in spite of the long- 
range interaction. More recently, the integrability of the class of quantum Hamiltonians 
has been shown, and algebraic structures for the underlying symmetry have been clarified 
systematically [26-37]. 

The asymptotic Bethe-ansatz solution (ABA) provides us with a systematic way to con- 
struct the energy spectrum [p] ,3]|, which has been known to give the exact solution to the 
continuum boson model and also to the SU(2) spin chain \TT)\. More recently, the ABA 
solution has been systematically generalized to multicomponent quantum systems with 1/r 2 
interaction, such as the SU(^) spin chain, the OSp(z/,l) supersymmetric t-J model, and a 
new hierarchy of the models related to the FQHE |12)| , |14)|Jl8)| , |21)|| . Furthermore, it has 
been found that a family of confined models with harmonic potential can be solved by the 
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renormalized harmonic oscillator (RHO) method [38~]|, which is a variant of the ABA. 

In this paper we wish to give a brief review of our recent studies on the critical properties 
of quantum 1/r 2 models based on the ABA solution and conformal field theory (CFT). This 
paper is organized as follows. In the next section we outline the ABA method by taking 
Sutherland's boson model as a simple example 0,0, and then discuss its critical properties 
based on CFT . In §3 the ABA solution is further extended to multicomponent systems 
with more complicated internal symmetry, the SU(V) Haldane-Shastry spin chain [|14)|] , and 
then in §4 the effects of hole-doping are discussed by the OSp(z/, 1) supersymmetric t-J 



model |12)| , |14)| which was introduced first by Kuramoto and Yokoyama for the OSp(2,l) 



case 



10] | . We further introduce in §5 a new family of 1/r 2 models which are closely related 



to a certain hierarchy of the FQHE ||18)| , [21)|| . In §6, we propose a new approach based on 
the RHO method ||38)|| in order to systematically construct the energy spectrum for the 
models with harmonic confinement |38)H42)|| . We then prove, by explicitly constructing the 
eigenfunctions for the grond state [|40)|| as well as the excited states |[41)| |, that the RHO 
indeed provides the exact spectrum for the above systems. Some applications to mesoscopic 
systems are also mentioned. 



§2. Critical Properties of Continuum Boson Model 



There are several variants of the integrable quantum models with 1/r 2 interaction. In 
this paper we will systematically investigate the systems with periodic boundary conditions 
and also the confined systems with harmonic potential. We note that quantum 1/r 2 models 



in continuum space was introduced by Calogero many years ago [Tj, and then have been 
actively studied by Sutherland p]| , p]| . In this section, we wish to mention characteristic 
properties of the quantum 1/r 2 systems by taking the boson model with periodic boundary 
conditions, and outline how to apply the ABA method. 

Let us introduce N interacting bosons in a ID chain of circumference L, 
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^=-eS+e^-^)> (i) 

with 1/r 2 interaction [7(r, — rj). In order to consider the model with periodic boundary 
conditions, the following form of the interaction should be used 



U(r) = V £ (r + nL)" 2 = Y£- sin~ 2 (^) = Vd~\ (2) 

n=-oo ^ 

which leads to V/r 2 in the limit of L — > oo. Here dij corresponds to the chord distance for the 
ring. It is known that the ground state for the model is given by the Jastrow wavefunction, 
i.e. the product of two-body functions 

^rilsinfc^ (3) 



with the Jastrow parameter, 



A= [VI + 2V + l]/2. (4) 



We will be concerned with the repulsive case V > hereafter. It should be noted that all the 
integrable 1/r 2 models take the Jastrow wavefunction as the ground state wavefunction. The 
two-body nature of the ground state is related to the integrability of the system, and implies 
that the two-body scattering may be essential for the many-body scattering. Based on 
this observation, Sutherland proposed the ABA method to exactly construct the excitation 
spectrum and full thermodynamics of the above boson system |2]] . 



2.1. Asymptotic Bethe-ansatz 

An asymptotic Bethe-ansatz (ABA) solution provides an elegant method to construct 
the excitation spectrum for the class of the above Hamiltonians Q. The essence of the 
idea is that although ordinary Bethe-ansatz (BA) methods are not applicable to the above 
systems with long-range interactions, the eigenfunctions can be written down consistently 
like those for the BA, 

V> = E A (Q; p ) ex p[* J2 k pj x Qj]i ( 5 ) 



in the asymptotic region, xq\ <C xqi <C • ■ • <C xqn, where Q (P) expresses one of N\ per- 
mutations for the coordinate (momentum) configurations. This form of the wavefunction 
implies that the many-body S-matrix in the asymptotic region can be factorized into two- 
body matrices consistently. If this is the case, one can diagonalize the many body scattering 
based on the factorized S'-matrix, and then obtain the spectrum of the system. The factor- 
ization of the S-matrix may be somehow suspected from the two-body nature of the Jastrow 
wavefunction for the ground state. At first glance the ABA solution seems to work only for 
low-energy excitations of continuum systems with the low density of particles. Remarkably 
enough, however, it turns out that the ABA solution exactly reproduces the whole energy 
spectrum for boson systems of any density j|2j |. Furthermore, this method has been found 



to be also applicable for the lattice models with high density of particles [|17) 



The two-body scattering in the Sutherland model yields the S'-matrix, Sy = 
— exp[—i9(ki — kj)] with the phase shift function 9(k) = n(\ — l)sgn(fc). Imposing pe- 
riodic boundary conditions, one can now diagonalize the many-body S'-matrix using the 



factorized S matrices, and then deduce the ABA equation for the rapidity kj [[2 



JV 



kjL = 2-Klj + (A - 1) ®( k j - h) (6) 

i=i 

with $(&) = irsgn(kj — fcj), where Ij is the quantum number which satisfies the selection rule 
Ij = (N + l)/2 mod 1. The interaction effects are now incorporated into kj via the phase 
shift function, and the energy is expressed simply in the noninteracting form, E = J2k^. 
Sutherland studied full thermodynamics based on the above ABA equations PK| . More 



recently, the critical behavior of correlation functions have been clarified ||11)|| with the help 
of CFT P1H46] 



2.2. Conformal properties 

By combining the ABA solution with the finite-size scaling in CFT, we now study critical 
properties of the model, and then evaluate the critical exponents of correlation functions. We 



briefly review the results obtained in ref. [IT)]. Let us start with low-temperature properties 



of the free energy [p)|, 



F(T) ~ F(T = 0) - ^ , (7) 
where the velocity of elementary excitations is t> = 27r\n. According to the known formula 



for finite-size scaling in CFT [ |44)|| , we can read the central charge of the underlying Virasoro 
algebra as c = 1 ||1 1 )|[ . Therefore, the critical behavior of the present model is expected 



to be described by c = 1 Gaussian CFT. To confirm this, we next calculate the finite-size 
corrections to excited states. Elementary excitations can be specified by the deviation of 
quantum numbers from the ground-state distribution. If we take Ij = (2j — N — m— l)/2 + d, 
the quantum number m labels the excitation which changes the particle number, whereas 
d denotes the excitation which carries the large momentum 2kpd with kp = itn. We then 



classify the excitation as [11 



27TU , , . 

AE ~ — — x(m; d; n ± ) , (8) 
Ij 

where the chemical potential term has been omitted. According to the finite size scaling in 
CFT |[45)|| , we can deduce the scaling dimension, 



x(m; d^n^) = ^m 2 + —d 2 + n + + n , (9) 

where particle-hole excitations are denoted by non- negative integers n ± . The selection rule 
for the quantum number reads: <i=integer for bosons. The momentum carried by the above 
excitation is 

2n 

P = 2irk F d +—[md + n + - n~] . (10) 

Ij 

From (9) and (10) we can determine conformal weights A 1 * 1 for the holomorphic (antiholo- 
morphic) piece, which characterize the operator content of the underlying Virasoro algebra 



in CFT |[43)|| . Recalling that CFT predicts the energy and the momentum to be universally 
related to conformal weights as x = A + + A~ and P = (2tt/L)(A + — A - ) for the 1/L sector 



1 45)|| , we obtain conformal weights ||1 1 )|| . 



A±{m-d-n±) = \{^±dR) 2 + n ± (11) 

with R = l/v 7 !. We now confirm that the above expression for conformal weights is inherent 
in c = 1 Gaussian CFT realized by free bosons with periodicity R, in which non-negative 
integers feature the conformal tower [46][. From the above analysis ||1 1 )| 



elude that the present model with 1/r 2 interaction is a typical example of Luttinger liquids 
47)H49)|| , i.e. its critical behavior is controlled by c = 1 CFT. 



2.3. Correlation exponents 

It is now easy to determine the critical exponents of correlation functions. For instance, 
let us begin with long-distance behavior of the density correlation function, 

< n(r)n(0) >~ const. + a r~ 2 + a 2 r~ a cos2fcpr , (12) 

where n(r) is the density operator. By taking (m, d, n ) = (0, 1, 0) as the quantum numbers, 
which carry the 2k p momentum transfer, we obtain the critical exponent (a = 2x), 

a = 2/A (13) 



for 2k p oscillation piece ||11)|| . Note that there is no logarithmic correction to the correlation 



function. On the other hand, the field correlator for bosons 



< 0l(r)0 6 (O) >~ r~ A , (14) 

has the leading non-oscillation term. By taking the set of quantum numbers (m, d, n^) = 
(1,0,0) for the primary field, we obtain the corresponding critical exponent ||1 1 )|| , 



Pb = A/2. (15) 

Fourier transform of this correlator yields the momentum distribution function around the 
origin, nf,(k) ~ \k\ 9b , with the corresponding critical exponent = — l = A/2 — 1. 

As the interaction strength V increases, the critical exponents vary continuously, charac- 
terizing the U(l) CFT critical line. For the special values such as V = and 4, the present 



exponents ||1 1 )|| agree with those obtained by Sutherland using the random matrix theory 
|2)| . We note that the system exhibits interesting properties at V = 4, at which the period- 
icity of bosons becomes R = 1/ y/2. This implies that symmetry of the model is enhanced 
to SU(2), and is described by the level-1 SU(2) Wess-Zumino-Witten model ||46) |. Therefore 



the effective theory for this special point is equivalent to that for the SU(2) Haldane-Shastry 
spin chain Q.Q. 

The above critical exponents satisfy the universal scaling relations for the Luttinger 



liquid, (3b = 1/a \17)\. It is instructive to note that in contrast to the known models with 



short-range interaction, the above critical exponents do not depend on the density of particles 



8),IT)|, but only on the interaction strength V. This peculiar property is characteristic of 
1/r 2 quantum models. We will encounter similar examples in the following chapters. 

We have been concerned with the boson case so far. The CFT analysis for the fermion 



case can be performed similarly |[L1) |. For fermions, the selection rule should read d 



m/2 mod 1, reflecting antisymmetry nature of the wavefunction. Hence the field correlator 
for fermions has the leading kp oscillation term and the corresponding exponent is given 
as Pf = [A + l/A]/2 by choosing the quantum numbers as {m.d.n^) = (1,1/2,0). This 
leads to the critical exponent of the momentum distribution around the Fermi point, Of = 
[A + 1/A-1J/2. 

2.4- Haldane-Shastry spin chain 

The quantum 1/r 2 model was extended independently by Haldane and Shastry to the 
lattice case, namely the antiferromagnetic S = 1/2 spin chain with 1/r 2 exchange interaction 



9)}. The Hamiltonian for a periodic ring with L sites reads, 



n = E Joists? + sfs] + i P (p - i)s;sfl , (ie) 

i<j 

where is chosen to satisfy periodic boundary conditions, J^- = Jd(xi — Xj)~ 2 with the 
chord distance d(x) = {L/tt) sin(7rx/L) (the i— th site is denoted by Xj). Here p represents 
the anisotropy of the model which is assumed to be p > 2 in what follows. The isotropic 
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model (p = 2) has been quite well investigated by Haldane and Shastry |8~)1,|9~)11 . Particularly 
in this case the ground-state wavefunction is given by the completely-projected Gutzwiller 



wavefunction at half filling [pj.9 



= exp(-ivr ^ x { ) d{xi - Xj) 2 , (17) 

i i<j 

in terms of coordinates {x{\ for down spins. As for the anisotropic case, the ground-state 
was also obtained in the Jastrow form for a positive even integer p Here we briefly 



summarize CFT analysis of the anisotropic Haldane-Shastry model fllTj 

It is known that the ABA solution is applicable to the above model although the asymp- 
totic region is not realized for a lattice system with high densities of particles |[L7)| . The 
two-body ^-matrix for the model is obtained as SV,- = — exp[— i8(ki — kj)} with the phase 
shift function 8(k) = (p — l)7rsgn(&;). Hence the ABA equation is essentially the same as the 
continuum boson case 117) 1, 



N 

kjL = 2txI 3 + (p - 1) £ Hkj - k t ), (18) 

i=i 

where N is the number of down spins. In contrast to the continuum model, however, the 
available range of kj is restricted to [— 7r, 7r], reflecting the periodicity of the lattice. One of 
the crucial consequences due to this restriction is that there exists the lower bound for the 
magnetization in the liquid phase, s z = 1/2 — 1/p. Therefore the massless phase is realized 
for magnetic fields, H c i < H < H c2 , where the upper critical field is given by 

F c2 =^[p(p-1) + 1], (19) 

at which the system is fully polarized, and the lower critical field is 

H cl = ^-pip - 1) - 1], (20) 

at which the magnetization takes its minimum value s z = 1/2 — 1/p in the liquid phase 
1 11)| . The ABA solution is applicable only for this range of magnetic fields. We would 



like to mention that the supermultiplet structure discovered by Haldane |17)| , which is now 



classified by the Yangian algebra ||36)|| , made it possible to construct the full thermodynamics 
for the isotropic case (jp — 2). 

Similarly to the Sutherland model, the low temperature free energy in the liquid phase 
takes the form, F(T) ~ F(T = 0) - 7rT 2 /(6v) with the spin velocity v = irp(l - 2s 2 )/(4J). 
This implies that the spin-liquid behavior of the present model is described by CFT with 
c = 1. Also, the excitation spectrum turns out to be classified in the same form as in (9). 
Hence, we can conclude that the spin-liquid phase for H c i < H < H c2 is classified as that 
of Luttinger liquids ||TTj 



Following the method outlined for boson case, it is straightforward to deduce correlation 
exponents. For instance, the asymptotic form of the spin correlation function is written as 

< S z (x)S z (0) >~ c x~ 2 + c 2 x~ a cos2k F x . (21) 

with the exponent a = 2/ p. On the other hand, the transverse spin correlation has the 
leading non-oscillation term 

< S + (x)S-(0) >~ x-P . (22) 

with the exponent (3 = p/2. These exponents indeed satisfy the universal scaling relation for 
the Luttinger liquid, a — 1/(3, as is the case for the continuum case. Note that conformal 
properties for the isotropic case (p = 2) has already been discussed by Haldane in detail 

eza- 



We have been concerned so far with the liquid phase with s z > 1/2 — 1/p under applied 
magnetic fields H c \ < H < H C 2. We then ask what will happen for H < H c \. Unfortunately 
the ABA solution is not efficient to answer this question. We can say that the system shows 
a singularity at H = H c i, but it is not clear whether the ordered phase is realized below H cl 
or not. We close this section by mentioning key points to solve the above question. We first 
point out that when the system approaches the boundary H c \ from the liquid phase, the 
spin correlation with the period p is enhanced at H = H c \ [fLlj ] . So, it may be interesting 



to check what kind of phase would be stabilized just below H&. It is also instructive to 
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ask what happens for the magnetization when H decreases further below H c \. According 
to the results in the Ising limit [ B] | , it may be possible that magnetization would show stair 
structures as a function of H. These interesting issues are to be clarified in the future study. 



§3. SU(i/) Spin Chain 

Now we wish to extend the quantum 1/r 2 models to multicomponent cases. Such a 
generalization was first made by Kuramoto and Yokoyama using the super symmetric t- 
J model, who found the Gutzwiller state as the ground state and discussed low-energy 
excitations ||10)|| . Subsequently, this model was solved by the ABA method |12) |, which has 



been later proven to give the exact spectrum of the model ||13)|| . Further generalization to the 
SU(z/) spin chain and to the multicomponent t- J model with OSp(z/, 1) supersymmetry has 
been done independently in refs. ||14)|J15)|| . Before discussing the supersymmetric t- J model 
we first study a SU(z/) generalization of the Haldane-Shastry model {y > 2), and then discuss 
the effects of hole-doping using the supersymmetric t-J model in the next section. 

We introduce the Hamiltonian of the SU(z/) spin chain with 1/r 2 interaction by the 
following general form [ |1 1)|| . 



n= l - J2 {-ifWJijXfxf*, (23) 

i<j,a,[3 

where the Hubbard operator Xf = \ia >< i(3\ interchanges states at 2-th site from (3 to 
a, and the exchange coupling Jy is given by inverse- square interaction (16). Here we have 
assumed that the system consists of v components of spins (or colors) with a, (3 — 1, 2, ■ • • , v, 
and the fermion numbers take F{(3) = (0, 0, • • • , 0) for (3 = (1, 2, • • • , v) in case of the SU(z/) 
spin chain. Hence all the particles obey the same statistics. This Hamiltonian, which is 
a SU(z/) generalization of the Haldane-Shastry model ||14)|J15)|| , is indeed invariant under 



the global SU(z/) transformation. If we add another state (a, = v + 1) with the different 
fermion number F(z/+1) = 1, this Hamiltonian describes the multicomponent t-Jmodel with 
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OSp(V, 1) supersymmetry which is realized by doping holes into the SU(z/) Haldane-Shastry 
model. 



3.1. Ground-state wave function 

We start by writing down the completely projected SU(z/) Gutzwiller wavefunction as 



the exact ground state for the above SU(i/) spin chain |14)| , |15)|| 



i^>=^ ) nn4: )t io>, (24) 

a=l k a 

where is the creation operator of electrons with spin a (1 < a < v) and momentum 
k a . We have assumed here that Pq selects configurations for which every cite is occupied 
by only a single electron with spin a. The projection Pq^ is easily done by taking the 
reference state \F > full of particles with z/-th spin [50j . In this representation, we rewrite 
the Gutzwiller state as, 

i^>=E^(K (Q) })n^ } v>, (25) 

(a,i) a,i 

where b) creation operator for spin particles (1 < a < v — 1). The Gutzwiller 

wavefunction is given in terms of coordinates of spin particles (x\ a ^), 

^G{{xt ) }) = zM^Y.* ( ? ) )T[d{x < f ) -xff I] d(x\ a) -xf), (26) 

where d(x) = (L/tt) sin(7rx/L) is the chord distance and L is the number of lattice sites. 
Following techniques developed for the SU(2) spin chain |8)1 , |9]| and the OSp(2,l) super- 



symmetric model IP- 0)|| , it is straightforward to show that the above SU(z/) Gutzwiller state 
without holes gives the ground state for the SU(y) spin chain | 14)| , |15)|| . The corresponding 
ground-state energy is thus computed as 

which will be shown to coincide exactly with the result deduced from the ABA solution. 
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3.2. ABA solution 

We now derive the ABA solution to the SU(z/) spin chain, following the calculation 
outlined in ref. [ 14)|| . Let us begin with the two-body scattering. In the asymptotic region 



in the coordinate space, the two-body scattering matrix for the above model is known to 
take the simple form, 

Sff = lim h ~ kj 7 + i7]Pa \ (28) 
J v-*o ki — kj — it] 

in terms the permutation operator P a p which interchanges the coordinates x a and xp. In- 
ternal SU(z/) symmetry is now simply taken into account via the operator P a p. It is quite 
remarkable that this form of the S'-matrix is essentially the same as that for noninteracting 
electrons. It should be noted, however, that this model actually describes a non-trivial spin 
system with long-range interaction. We note that Sfj satisfies the Yang-Baxter factorization 
equation, Oj k o ik — o ik jj k . 

Consider now the many-body scattering among particles with v — 1 different spins by 
taking the i^-th species as the background. As mentioned above all the particles obey the 
same statistics for the present SU(z/) spin chain. Imposing periodic boundary conditions for 
the ring system with L sites, we have to diagonalize the scattering problem 

<" :j '* S {j+1)j S {mj ---S {j (29) 
where SL- = S l A ||51)|| . Introducing v — 1 kinds of rapidities kj for a — 1, 2, - • • , v — 1, one 



can solve this problem by the nested BA method ||51) , 52)| |. We finally arrive at the nested 



ABA equations for the rapidities |14 



kfL = 2nl^ + £ S(*g> - Af ) + J2 - kj\ (30) 

m i 

E - k?) + 2*1$ = E $(4 2) - kf) + E *(4 2) - kf\ (si) 



e *{* { r x) - k [ r 1] ) + 2vr/r i) = e - ^~ 2) ), (32) 
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with <&(k) = 7rsgn(fc), where 1^ classifies v — 1 kinds of spin excitations, 



jj Q) = \[M V _ X + M U + M v+l ) mod 1, (33) 
with Mq = M v = 0. Here have introduced the quantity 

M a = J2Np, (34) 

f3=a 

in terms of the number of particles N a with spin a = 1, 2, • • • , v. Henceforth v — 1 kinds 
of spin excitations are referred to as spinons for simplicity. The energy is given in a non- 
interacting form, 

Mi i 

E = T,-M ] ) 2 -- 2 } + a (35) 



with the energy shift 



a = g(l-^)(L-2M 1 ) + ^(l--^). (36) 



3. 3. Bulk properties 

We now introduce the density function p a (k) for the spin rapidity k^ in the thermo- 



dynamic limit ||14)|| . Note that the density function p a {k) has constant values in the region 



[R a : B a+1 < \k\ < B a ] with the condition Bi > B 2 >,••-,> -B^-i and B v = 0. These 
density functions are easily evaluated as, 

= < I 37 ) 
I 0, /3 < a 

for the region R@. We note that p a has a discontinuity at every boundary B a . The density 
of particles with spin a is expressed in terms of p a (k) as 

/B a rB a+ i 
p a (k)dk — / p a+ i(k)dk (38) 

for a = 1, 2, • • • , v — 1 and that for the u-th spin is determined by the sum rule: YX%=\ n a — 1- 
The cut-off parameter 5 a is controlled by the applied magnetic field through the relation, 

14 



B a = [tt 2 - 2a{a + 1)H} 1/2 , a = 1, 2, ■ ■ ■ , v - 1, (39) 

if we assume that v kinds of states have magnetic moment S, S — 1, • ■ • , — S with v = 2S + 1. 
Note that B a = n corresponds to the SU(z/) singlet state while B a = to the fully polarized 
state. 



It is now straightforward to compute bulk static quantities [ [14)11 • For example, the energy 
for SU(i/) singlet is given as 

which agrees with the result (27) obtained by the SU(z/) Gutzwiller wavefunction. Also, the 
magnetization is evaluated as a function of magnetic fields, 

s z = V —±-^ X> 2 - 2a(a + 1)B\W0(H<?> - H), (41) 

2 27T , 

where 6{x) is a step function and we introduced the critical field by Hj^ a ' = n 2 /[2a(a + 1)]. 
This yields the uniform spin susceptibility at low fields, 

Xs = tL E«(« + 1 )- ( 42 ) 

It is seen from (41) that the contribution to the magnetization from a-th spinons is saturated 
at H = H( a \ and the effective degrees of freedom for spinons are reduced successively as the 
magnetic field is further increased. Therefore spin-liquid phases realized in magnetic fields 
are classified into v — 1 sectors for which each phase boundary is determined by H^ a \ It is 
also easy to evaluate the velocity of each spinon excitation, which is given by the formula, 

v a = ^[n 2 -2a(a + l)H} 1 / 2 . (43) 

This velocity determines the specific- heat coefficient, 

C/T = IJ2-6(H^-H). (44) 

6 a=l V a 



According to the finite-size scaling in CFT ||44)[| , the above expression for C/T implies that 
v — 1 kinds of spinon excitations are described by independent c = 1 CFT. This point will 
be discussed in detail below. 
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3. 4- Conformal properties 

We now classify low-energy spinon excitations in order to observe conformal properties 
[ ■ Introduce first two vectors in and d out of quantum numbers associated with v — 1 



kinds of spinon excitations. A quantum number m a (a — 1, 2, • • • , u — 1) labels the change of 
the number of a-th spinons, i.e. m a = AM a . On the other hand, d a denotes an excitation 
which carries a momentum 2nM a /N . Hence the total momentum is expressed as, 

v-l 

AK = 2(tt - ak F )d a , (45) 

a=l 

where kp = tt/u for H = 0, and 1/L corrections to the momentum transfer have been 
neglected. According to boson nature of spinon excitations, the quantum numbers d a satisfy 
the selection rule 

d a = ^K-i + m a+1 ) mod 1, (46) 

for 1 < a < v — 1 with the condition m = vn v = 0, which directly follows from eq. (33). 
Low-energy excitations are now classified in the matrix formula, 

2,71 

AE = —Y, v a x a , (47) 

L a=l 

where x a is the scaling dimension for each spinon with spin a, 

x a = \{Z- l rh)l + {T}d)l + ni + n", (48) 

and labels particle-hole type excitations for the a-th spinons. The [y — 1) x {y — 1) 
matrix Z, which is referred to as the dressed charge matrix [53]], is given as [|14)| 



' a/3 



aW + l)}-"\ a</3<v-l, 

(49) 

0, a > (3. 



One can see that the above scaling dimension has the typical form inherent in c = 1 CFT 
46)|| and that non-negative integers form conformal towers which characterize the rep- 



resentation of Virasoro algebra ||53)|| 
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Let us now consider the total scaling dimension [|14) 



x 



E 

a=l 



X r 



(50) 



which determines the critical exponents r] = 2x for correlation functions. It turns out that 
x is expressed in terms of the Cartan matrix C for SU(^) Lie algebra, 



1 - "~ 1 

x = -rtfCm + iC~ 1 d+ J2( n t + n a), 

^ a=l 



(51) 



where 



/ 2 -1 
-1 2 -1 



(52) 



V "12/ 

A remarkable point is that this expression for x holds even in magnetic fields although there 
does not exist SU(z/) symmetry and the velocities of spinons are strongly modified. These 
characteristic properties are in contrast to those for the SU(z/) antiferromagnetic Heisenberg 
chain with nearest neighbor interaction for which the scaling dimension changes continuously 



and the Cartan matrix can specify excitations only at H = ||54)|| . Therefore the correlation 
exponents are not affected by magnetic fields in the present model, and are determined by 
those of c = v - 1 CFT. 



3.5. Correlation exponents 

Following standard methods in CFT [ |43)H45 )|| , we deduce all the critical exponents of 



various correlation functions |14)|. We apply the assignment of quantum numbers to a given 



correlation function, which was used for the Hubbard model [481] and subsequently applied 



to the t-J model with nearest-neighbor hopping [ [49 )|| . For example, let us consider the 
long-distance behavior of the spin correlation functions characterized by the exponent f3j, 



< S z (r)S z (0) >~ A J cos(2jk F r)r~ 



•Pi 



(53) 
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This correlation function conserves the number of particles, so that we set m = (0, • • • , 0). 
There remain several choices for the quantum numbers d. A choice of d a = 5 a j (1 < j < v—1) 
results in the 2jkp oscillation piece of the spin correlation functions. We thus obtain spin- 



correlation exponents for the 2jkp oscillation part [14 



ft = ^4 1<J<»-1. (54) 
For the SU(2) case, this result agrees with that of the Gutzwiller wavefunction at half filling 



59)|| . We stress here again that the above critical exponents do not depend on magnetic 
fields. 



§4. OSp(z/, 1) Supersymmetric t- J Model 

Now we wish to observe what is modified when we dope holes into the SU(z/) Haldane- 
Shastry model. This problem seems interesting since highly correlated electron systems have 
attracted particular attention recently. Such an attempt was firstly done for the OSp(2,l) 
supersymmetric model ||10)|| , and subsequently the model was extended to more generic 
supersymmetric model with OSp(V,i) symmetry [|14)| , |15)|| . We give a brief review of ref. 



141 here. 



4 ■ 1 ■ Symmetry properties 

We first note that if the fermion numbers are taken as F{(3) = (0, • • • , 0, 1) for (3 = 
(1, 2, • - • , u+ 1) in the Hamiltonian (23), this model is equivalent to the multicomponent t-J 
model PIPl , 



« = -EE ^Ua + Ey[ E " (1 " "*)(! - %)]> (55) 

i<j a=l i<j a,/3<u 

with the so-called supersymmetric condition tij = Jij, where ni is the electron number at 
the i-th site. Here the double occupation of every site is strictly prohibited. We have 



used the fact that the term like x^ u+l ^ xjj v+1 ' a with a = 1, 2, • • • , v corresponds to the 
hopping term of electrons with spin a. For a special case with three-component particles 
with F = (0, 0, 1), the model reduces to the OSp(2,l) supersymmetric t-J model (apart from 
the chemical potential term) 110)11 , 



n = - Yl '■j'Wm + Yl M s i ■ s i - j n i n jl ( 56 ) 

i<j,a i<j 

which is obtained by doping holes into the SU(2) Haldane-Shastry model. 

Here, we briefly mention supersymmetry properties of the above t-J model. For the 
graded model (23) or (55) with bosons and fermions, it is known that the Hubbard operator 
should satisfy the following commutation relation |p5)|,|56)|| . 



[Xf, X™\ ± = S^xr ± S av Xf), (57) 

where the anti-commutator (+) should be used only for the fermion operators. This algebra 
is called as the doubly-graded Lie superalgebra with Osp(z/, 1) supersymmetry. Therefore 
the model (55) is invariant under the global OSp(z/, 1) transformation |p~5) , |56)|| . Thanks to 



the supersymmetry of the Hamiltonian, we can successfully apply the ABA method to the 



above t-J model ||57) . 



4-2. Ground-state wavef unction 

Before proceeding with the ABA calculation, we mention the ground-state wavefunction. 
The ground state of the OSp(f, 1) model is given by the SU(z/) Gutzwiller state with holes, 
as was demonstrated for OSp(2,l) case ||10)|| . In this case the Gutzwiller projection operator 



in (24) projects out configurations with more than one electron on each lattice site, which 
implies that there can be cites without any electrons (assigned as holes). Taking the fully 
polarized state \F > as the reference state, the Gutzwiller state with holes is written as 



\* G >= e *Ga*PMss}mb<r»i[h]\F>, m 
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where ty is a creation operator for spin particles and hi = is that for doped holes. 
The Gutzwiller wavefunction is now expressed in terms of coordinates of spin particles (x\ a ^) 
and holes (s m ) as, 



* G (K (a) },{^}) = expHvr(E^ ) + E^)] II d{^-^ 



a,i<j 



x n d(x\ a) - xf) n d(4 a> - *») n ^ - «„) (59) 

a</3,i,j a,i,m m<n 

One can show easily that this gives the ground-sate wavefunction of the OSp(z/, 1) super- 
symmetric t-J model, the detail for which can be found in refs. |T0), 14] , |15)| . 



4.3. ABA solution 

We shortly outline how to deduce the ABA solution to the Osp(z/, 1) supersymmetric t-J 
model ||14)|| . A remarkable point is that the two-body S- matrix for this model is given by the 
same formula (28) as in the SU(z/) spin chain, and symmetry properties can be taken into 
account solely via the permutation operator P a/3 . Hence, using the ^-matrix (28) for the 
two-body scattering, it is now straightforward to deduce the ABA solution to the OSp(z/, 1) 
t-J model. Taking the reference state full of particles with the ^-th spin, consider now the 
scattering problem among v — 1 kinds of spinons and holons (holes). In order to diagonalize 
this problem, therefore, besides v — 1 kinds of spin rapidities (a — 1, 2, • • • , v — 1), it is 
necessary to introduce the charge rapidity kj ||52)|| . As a result the last line of the ABA 



equations in (32) should be modified by grading holes [I4~)|, 



k?L = 2nl^ + £ $(4 2) - fcf) + E - k\\ (60) 



E *(*£> " kf ] ) + 2vr4 2 ) = E *(4 2) - kf) + E *(*S> - ki% (61) 



27T/M = E *(kM - kf~\ s = 1, 2, ■ • • , M„, (62) 
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where I^"' is a quantum number which labels the degrees of freedom for doped holes. Ac- 
cording to the antisymmetric nature of electrons, the selection rule for this quantum number 

reads, 

= -M„ mod 1. (63) 



4-4- Bulk properties 

We calculate bulk quantities following techniques outlined for the spin chain. It is re- 
markable that the shape of the the calculated magnetization curve does not depend on 
electron concentrations, which is given by the same formula as (41), but only the critical 
field 

Hie) = *Mn~2) 

2a(a + l) 1 1 

is modified by hole-doping (n: electron density) ||14)|| . Also, other bulk quantities are easily 
calculated. By the second derivative of the ground state energy with respect to n, we obtain 
the charge susceptibility (compressibility) as a function of electron concentrations, 

2u . . 

Xc = -^f. r, 65 

7i 2 (l — n) 



which shows a divergent behavior near insulating phase n = 1 [rO],I4),15j|. The coeffi 



cient of the T-linear heat capacity is given in terms of the velocities of spinons and holons 



7T r l ^ 1 



7 = -[- + £ -£(#(«) -#)], (66) 



3 v c ^ v a 



where the velocity of spinon excitations, v a , takes the same formula as in (43) while that 
for holon excitations depends on the electron concentration, 

v c = n(l-n)/2. (67) 
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The above expression for the specific heat implies that the critical behavior of the present 
model is described by c = 1 CFT. We can see that there are the Luttinger-liquid relation 
among bulk quantities, vrx c f c = v. 



4 ■ 5. Conformal properties 

We now classify low-energy spin and charge excitations in order to study conformal 
properties ||14)|| . Let us introduce //-component vectors fa and d for quantum numbers which 



specify excitations. We assume the v-tYt component of vectors to be related to the holon 
degrees of freedom. According to the antisymmetry properties of electron wavefunction, the 
selection rule for quantum numbers reads 

da = ^(m ft -i + m a+ i) mod 1, (68) 

for spinon excitations (a = 1, 2, • • • , v — 1), whereas that for holon excitations is 

d y = -(m u _i + m u ) mod 1. (69) 

Note that quantum numbers d a carry the large momentum transfer 2(ix — akp)d a for 1 < 
a < v with Fermi momentum kp = Tin/ v. The excitation spectrum is now classified as 

AE = —Y / v a x a , (70) 

L a=l 

from which we can read the scaling dimension x a . The resulting scaling dimension is ex- 
pressed in the same formula as (48) by extending the dressed charge matrix ||53)|| to the v x v 
matrix Z, 

' aW + \)\-y\ a<(3<u-l 7 

a/y/u, (3 = v, 1 < a < v 



and Z a p = otherwise [|14)|| . Matrix elements of u-th row and the u-th. column are related 
to the holon degrees of freedom. From the expression (70), one can see that that holon 
excitations as well as spinon excitations are described by independent c = 1 CFT. Namely 
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c 



(72) 



this liquid is classified as the Luttinger liquid. The total scaling dimension x = YZt=i x a is 
reduced to the simple expression (50) where the v x v matrix C in this case is 

/ 2 -1 \ 
-1 '•• -1 

'•• 2 '• 

V -i i 7 

which is nothing but the Cartan matrix for the OSp(i/, 1) Lie superalgebra |14)| . We wish to 
emphasize here again that the above scaling dimension depends neither on magnetic fields 
nor on electron concentrations, characterizing the universality class of quantum 1/r 2 models. 



4-6. Correlation exponents 

Let us now evaluate critical exponents for various correlation functions with the use of 
CFT. First we compute the critical exponents of the spin correlation function. By choos- 
ing the quantum numbers d a = 5 a j (1 < j < v — 1) for the spin correlation function 
48lPlj57l|58| , we obtain the spin correlation exponents for the 2jkp oscillation part, 



Pj = 2j, l<j<u-l. 



(73) 



Remarkably enough, all the exponents are given by the canonical (integer) values. Also, 
critical exponents for the charge correlation functions are obtained similarly, which are 
given by the same formula as for the spin exponents, but the 2vkp oscillation piece appears 
with (3 V = 2v. It is instructive to note that the spin correlation exponents discontinuously 
change when the holes are doped into the spin chain. This can be checked by comparing 
the expressions (54) and (73). This property of the discontinuity is known to be common in 
the correlated electron systems close to the insulator. 

We now evaluate the correlation exponent rj of the field correlator of electrons, < 
c a( r ) c a(0) >~ cos(/ci?r)r _r? . A simple excitation relevant to this correlator is given by the set 
of quantum numbers (mi, m 2 , • • • , m v ) = (1, 1, • • • , 1) and (d\, c^, • • • , d„) = (1/2, 0, • • • , 0) 
|18)| , |49)| , |57]1 , |58)|| . The momentum carried by this excitation is 7r(l — kp). The resultant 
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critical exponent turns out to be 77 = 1 ||14)|| . Fourier transformation of this correlation 
functions gives the momentum distribution 



Tik = Uk F — const. \k — kp\ e sgn(k — /cp), 



(74) 



with the corresponding exponent 9 = rj — 1 = 0. By taking account the fact that there is no 
logarithmic correction in the present case, 6 = implies that there should be a discontinuity 
in the momentum distribution at the Fermi point, as firstly pointed out for the OSp(2, 1) 



model |[L0)|| . Recall again that all the critical exponents of correlation functions are given 
by the canonical values as well known for the Gutzwiller wavefunction [|59)| . At first glance 
these canonical exponents seem to be contradicted to Luttinger liquid theory ||17)|| , but we 
should keep in mind that this fixed point is indeed on the critical line of the Luttinger liquid 
(c = 1 Gaussian CFT). The present model is, therefore, classified as a specific example of 
the Luttinger liquid which has canonical exponents ||10)| . 

We have not been concerned here with thermodynamic properties at arbitrary temper- 
atures. In order to describe the full thermodynamics of the lattice models, it is crucial to 
study the degeneracy of excited states. Concerning this problem, the free-spinon picture of 
Haldane [|17)|| , Yangian symmetry ||36)|| and related methods ||13)| . |32)|| have been successfully 
used to construct correct thermodynamics. 



§5. Hierarchical Models Related to FQHE 



One of the most interesting aspects of the 1/r 2 quantum systems is the intimate rela- 
tionship ||17)H22)|| to the fractional quantum Hall effect (FQHE) [|60)H62)|| . We have indeed 
seen that the construction of eigenstates for the 1/r 2 models is quite analogous to that for 
the FQHE: the ground state is given by the Jastrow wavefunction, and the excited states 
are constructed by multiplying polynomials to the ground state wavefunction. Therefore, 
besides much interest in the integrability, these quantum models should exhibit interesting 
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phenomena related to the FQHE. We have recently proposed a novel hierarchy of the ID 
quantum models with 1/r 2 interaction ||18)|| , the construction of which is essentially same as 
that of a hierarchical FQHE with the filling fraction ||61) , |62)|| , 

1 



fu 



1 



(75) 



Pi 



Pr-- l 

Pv 

In particular the matrix deduced from the energy spectrum has been shown to coincide with 
the topological- order matrix which characterizes the internal structure of the FQHE state 
||63)H65)1 . In this section we briefly review the hierarchical models related to the FQHE 



181,21; 



5.1. Continuum Models and ABA solution 

Let us introduce a family of ^-component electron models (a = 1,2, ■ ■ ■ ,u) with 1/r 2 inter- 



^ = "2E^+EE d^V a ,(V aP + P°f), (76) 



action in the periodic ring of length L |]T8 

1^ d 2 

a<P i<j 

with chord distance d%j = (L/tt) sin[7r(xj — Xj)/L], where P°f is the spin exchange operator 
with spin indices a, (3 = 1,2, ■ ■ ■ v, and V a p is the interaction parameters dependent on spin 
indices. Starting from noninteracting SU(^) electrons (zero-th stage), let us specify the 
stage of the hierarchy in terms of the matrix formula of the interaction parameter. The first 
stage is defined by introducing interaction Ai uniformly to all the electrons, i.e. V a p = \\. 
The resultant model is the SU(^) Sutherland model for which only the charge excitation is 
affected by the interaction ||15)| , |27)| , ^8]| . The second family is introduced by turning on the 
uniform interaction V a p = A2 among electrons except for the a = 1 species. By iterating 
this procedure step by step, the v-th stage of the model is characterized by the interaction 
matrix ||18)|| , 
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/ Ai 



V 



Ai\ 



VAi 



+ 



Ai/ 



A 2 



\ 

A 2 



/O 



+ • ■ ■ + 



o \ 



o KJ 



(77) 



\0 A 2 •■■ x 2 j 

Note that the interaction A m acts on particles with spin indices a = m, m + 1, • • • , v, hence 
it is regarded as spin-dependent interaction. The construction of the above hierarchy is 
quite analogous to that for the FQHE ||61)| , |62)|| , in particular, to Jain' composite- fermion 



construction of the hierarchical FQHE ||62)|| . In the FQHE case, the interaction matrix 
introduced here serves as a topological order matrix which specifies the internal structure 
of the FQHE ||63)| - p5)|| . We shall see that the interaction matrix completely specifies the 
low energy excitations, and that the present model can describe essential properties for edge 
states of the hierarchical FQHE with filling fraction f v ||63)| . 



Let us now find the solution to the problem by the ABA method |18)| ] . We first consider 
the scattering of two particles by the interaction d^V a p(y a p + 1) As mentioned in the 
previous sections, this interaction yields the S-matrix of exp[— i$ a p{ki — kj)] with the phase 
shift function, 



<fiap(k) = V a p%sgn(k), 



(78) 



in the asymptotic region (|xj — Xj\ >> 1). The step-like form of the phase shift <fi a p{k) is 
characteristic of the 1/r 2 interaction. Taking into account the spin degrees of freedom, we 
then find the two-body S- matrix for the above hierarchical models as [FS]], 



Su = lim k \ k \ ^ 
e^o hi — — ie 



aj3 



-i^alliki-kj) 



(79) 



where the first factor arises from non-interacting SU(z/) electrons, while the second is due 
to the 1/r 2 interaction. The key feature that the S-matrix (79) is a simple superposition 
of two S'-matrices makes it possible to treat the scattering problem consistently by nested 
Bethe ansatz techniques. 

As mentioned in the previous section, the essence of the ABA method is that the many- 
body S'-matrix in this family can be decomposed into two-body matrices in spite of long- 
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range nature of interaction. For a simple case of SU(z/) symmetry, we have already checked 



that this solution gives the exact spectrum ||14)| , f21)|| . If this is the case in general, we can 
diagonalize the many body S- matrix by the standard nested Bethe ansatz. Consequently 



we arrive at the ABA equations for z/-kinds of rapidities k" 1 IT 



M 2 Mi 

kf>L = 2nl^ + £ $(4 2) - kf) + A a £ *(kW - ki\ (80) 

m=l 1=1 

M a M a + s 

{K + 1) E - k^) + 27T/W = E E ^ - kf + \ (si) 

1=1 s=±l j=l 

for 2 < a < v, where $(A;) = 7rsgn(A;) and 1^ is an integer or a half integer which classifies 
the charge and spin excitations. In the above equations we have introduced the quantity 
M a = Yl^= a N/3 where Np is the number of electrons with spin f3. The energy is written in 



the noninteracting form E = (1/2) ) 2 - 



5.2. Bulk properties 

Let us calculate bulk quantities ||18)|| . It is remarkable that all the bulk quantities are 
expressed solely by the parameter f v introduced in (75) if we assign the set of parameters 
Pi in terms of interaction parameters Aj as, 

Pi = Xi + 2- 5 lU (82) 

for a = 1, 2, • • • , v. For instance, the ground-state energy is evaluated as 

E a /L = ** n *IWl), (83) 

where n is the electron density. The second derivative of E g {n) yields the charge suscepti- 
bility (compressibility) in terms of f u , 

Xc = (fuM 2 n-\ (84) 

Roughly speaking, the quantity 1/ f u corresponds to the size of the exchange-correlation hole 
due to the 1/r 2 interaction: we can regard the effective volume of electrons to become 1/ f v 
times as large as that of free electrons. 
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The free energy at finite temperatures (T) is formulated as, 

/oo 
log[l + exp(-ei(Jfe))]dife 
-oo 

in terms of the dressed energy, 

ei(A;)/T = \k 2 - H - log[l + exp(-e 2 (fc))] + \ x log[l + exp(-ei(Ar)/T)], 



(85) 



(86) 



e a (k)/T= (-l) s (l + XJ q0 )log[l + exp(-e a+s (k))/T], 

s=-l,0,l 



(87) 



for a = 2, 3, • • • , u, with e v+ \ = oo. The low-temperature expansion of the free energy yields 
the coefficient of the T-linear specific heat as, 



C/T 



71 

3v 



(88) 



with the velocity v = 7rn/(2/„). Note that all the velocities for v kinds of elementary 
excitations have the same value v although there is not SU(z/) symmetry in the model. The 
Luttinger-liquid relation between the charge susceptibility and the velocity takes the form, 
ttXcV = /i//2, which can determine the critical behavior of charge excitations. 

5.3. Excitations 

Let us now turn to the excitation spectrum [18]]. Using the ABA equations, low-energy 
excitations are classified in the matrix formula, 

2vrw r l 



L L 4 

where the v x v matrix T is evaluated as, 

( Pi -1 
-1 Vi 



\ 



[-m'Tm + d'iT)- 1 ^, 



(89) 



(90) 



-1 PuJ 

Here the ^-component vector m is out of quantum numbers which classify the charge and 
spin excitations, where we assumed that mi and d\ label charge excitations. It is remarkable 
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that the above matrix T deduced from the energy spectrum is nothing but the topological- 



order matrix for the FQHE with filling fraction f v |64)|. One can see that there is v free 
parameters, Pi, in the above matrix. This implies that the critical behavior of the present 
model is described by the corrections of v independent CFT with the central charge c = 1, 
where Pi can make c = 1 critical lines. 

5.4- t-J model and lattice effects 

To see the relationship to the composite fermion theory in the FQHE more explicitly, let 



us observe what happens for the lattice case [ |18)|| . We introduce a hierarchical family of t-J 
models with 1/r 2 interaction, 



^=E^ 2 4c JQ + £ dgVpiVp + Pg 1 ), (91) 



with interaction parameters V a p defined in (77). Here configurations with more than one 
electron at each site are assumed to be prohibited. The first family of the hierarchy coincides 
with the SU(z/) t-J model introduced in | |15)|| . This class of lattice models can be solved by 



the ABA method, and the resulting nested ABA equations are given by the same formula 
as in (80) and (81). Hence, the bulk quantities are given by the same expressions as for the 
continuum case. In the lattice t-J model, however, we encounter two crucial constraints |[L8)|| : 
(a) parameter V a p should be an even positive integer, and (b) the rapidity kj should be in 
the region [—it, vr] as is the case for ordinary lattice models. According to (a), the parameter 
f v should be a fraction with the odd denominator, which demonstrates the analogy to FQHE 
explicitly. Furthermore, the constraint (b) brings about a remarkable property, i.e. a singular 
property at the electron density n = /„ above which the Luttinger-liquid state breaks down 



( n > fv) [ 18)|| . This means that a band-edge singularity for non-interacting SU(^) lattice 
electrons at the density n = v is modified into the singularity at the fractional filling n = f v 
in the presence of the 1/r 2 interaction. We note that this phenomenon is essentially the 
same as that for the composite fermion theory of the FQHE, i.e. the hierarchical FQHE 
with the filling n = f v can be realized by starting from the integer (noninteracting) QHE 
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with the filling v 



5.5. Chiral constraint and FQHE edge states 

In order to observe how remarkably the present model reproduces essential properties 
expected for edge states of the FQHE, let us think of what will happen if we consider only 
right (or left)-moving electrons in the ID system ||18)| , |21) |. This is referred to as chiral 



constraint which is essential for edge states of the FQHE in disk geometry ||63)|| . In order to 
deal with chiral constraint, we make use of a trick valid for the 1/r 2 systems. We first add 
two electrons at the left and right Fermi points, in order to suppress the current d which 
does not exist in the chiral model. This results in the energy increment both for the right 
and left branches. Dividing the energy increase into two parts and discarding the right-going 
piece, we obtain the excitation spectrum for left-going piece as | |18)| , |21)[ | , 

e = — m'Tm. (92) 
L 

This formula reproduces the spectrum described by the holomorphic piece of v independent 
c = 1 CFT. Hence we can determine the critical exponents for various correlation functions. 
For example, the critical exponent 9 a for the momentum distribution function 

njf } ~ const. + a sgn(A; -k F )\k - k F f a (93) 

is obtained as 

Q = X>. ( 94 ) 

All the other exponents can be similarly obtained. We stress here that all the critical 
exponents agree with those of effective field theory (chiral Luttinger liquids) for the edge 



states of FQHE with the fraction f u ||63)|| . In particular the key matrix T deduced from the 
excitation spectrum coincides exactly with the topological order matrix which characterizes 
the internal structure of the FQHE with the filling fraction f v ||64)|J65 j 



We have observed that the hierarchical models considered here describe characteristic 
properties for the edge states of the FQHE with filling fraction f u so remarkably. The 
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correspondence is not accidental, and we can indeed see a clear reason for it in the role 
played by the phase shift function A Q 7rsgn(fc). Recall that in the composite fermion theory 



for the FQHE ||62)||, flux quanta are attached to electrons in order to evolve the hierarchy 



starting from noninteracting electron systems [62)|- 65)|| . A crucial point is that attaching 
X a flux quanta in the FQHE corresponds to introducing the phase shift function X a irsgn(k) 
in the present ID system, which can be given by the 1/r 2 interaction | |18)| |. Hence, one can 
clearly see from this correspondence why the family of 1/r 2 models can describe characteristic 
properties of the hierarchical FQHE remarkably. 



5.6. Dual bases: holon-spinon and electron 

We would like to point out another instructive relationship to the FQHE ||21)|| . So far, 
we have classified the excitation spectrum in terms of the matrix T and the corresponding 
quantum numbers out of charge and spin excitations. This basis, which is referred to as 
holon-spinon basis, always shows up when we use the Bethe ansatz method. We should 
recall that there is an alternative basis, i.e. electron basis. Although the electron basis may 
not classify the excitation spectrum when various excitations have different velocities, it still 
describes the critical behavior of correlation functions correctly ||21)|| . Quantum numbers in 
the electron basis are obtained from those of holon-spinon basis via a linear transformation, 



with the matrix 



N = Um, J = (U 



(I -1 

1 "• 

U 



V 



t\-i 



d. 



(95) 



-1 
1 J 



(96) 



The matrix T is then transformed into the v x v symmetric matrix pi) 



T = I + V, 



(97) 
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where I is the v x v unit matrix. Remarkably, we encounter the matrix V out of the 
interaction parameters (see (77)). Hence the matrix V itself can classify the excitation 



spectrum in the electron basis ||21)|| . We should emphasize here that there exists the exactly 
same matrix (97) in the 2D FQHE which also characterizes the internal structure of the 
hierarchical state [|S4j| , |65j[| . The above two kinds of bases are called the symmetric basis and 
the hierarchical basis in the FQHE, both of which characterize the same topological order 
of the hierarchical FQHE [|64)| , |65)| . 



5.7. Example: SXJ(v) Sutherland model 

As a simple example, we consider the first family of the hierarchical models (76), i.e. the 
SU(z/) Sutherland model characterized by the parameters V a p = A ||15)| . |27)| . p8]|] . In this case 
the exact wavefunction has been obtained in the Jastrow form ||15)|| , 



^=ni4 a) -4?iVo 

l>m 

Here ipo is the wavefunction for SU(z/) free electrons, 



(98) 



Ifo = W* ] \ 2da II - z[ i ] ) 5aiam exp^vr sgnfa - a m )\. 



(99) 



where zffi = exp(27ra;^ / L) , L is the length of the periodic system and xffl are spatial coor- 
dinates of electrons with spin a (= 1, 2, • • • , i>). The current 2d a carried by a-spin electrons 



is assumed to take an integer value ||15)|| . This form of the Jastrow wavefunction clearly 
demonstrates the analogy to Jain's construction of the wavefunction for the FQHE with 



filling factor /„ = z//(z/m + 1) ||62)|| . Namely, starting from the noninteracting wavefunction, 
interacting electrons can be described by introducing the Jastrow factor. In the electron 
basis, the excitation is classified by the symmetric matrix [15],2T]1, 

/1 + A ••• 1 + A\ 



VI + A 



(100) 



1 + A/ 
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In order to observe the symmetry property, it is more convenient to make use of the holon- 



spinon representation, which is characterized by the following matrix \\21) 

/1 + A -1 \ 
-1 2 -1 



(101) 



V "12/ 
From this expression we can see U(l) symmetry for the charge sector (T n ) and SU(y) 
symmetry for the spin sector (T a! 3 for a > 2). Particularly the {y — 1) X (i> — 1) matrix out of 
T Q/ g for a, (3 > 2 is the SU(^) Cartan matrix which characterizes the Lie algebra with SU(^) 
symmetry. Hence, the critical behavior of the spin sector is characterized by level-1 SU(^) 
Kac- Moody algebra with the central charge c = v — 1, and the holon sector is governed by 
c = 1 CFT for which the scaling dimensions vary continuously according to A. 

It is now easy to obtain critical exponents for the correlation functions. For example, 
the critical exponent for the momentum distribution function is obtained as ||21)|| 



(102) 



which is in contrast to the chiral case, 



l c) = A. 



(103) 



Note that the critical exponent for the chiral case exactly coincides with the result for chiral 



Luttinger liquids of the edge states with filling fraction f u = vj{y\ + 1)) ||63)||. 



§6. Confined Models and Renormalized Harmonic Oscillators 



We have been concerned so far with quantum models with periodic boundary condi- 
tions. There is another class of the integrable 1/r 2 models with harmonic confinement 



1) 


21,29) 


34) 


.35) 


38) 


-41) 
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34)1 , there have also been several attempts to apply them to conductance oscillations in 
mesoscopic systems ||66)| . p7)| . |39]| . [40)|| . In this section, we propose a systematic construction 
of the energy spectrum for the class of 1/r 2 models with harmonic confinement ||38)|| . This 



approach is referred to as the renormalized-harmonic oscillator (RHO) solution, since the 
essence of the idea is that all the interaction effects are incorporated in terms of the renor- 
malized quantum numbers of oscillators. In this sense, this method is regarded as a variant 
of the ABA. We use this idea for a systematic construction of the spectrum for the confined 
1/r 2 models [|38)|| . By constructing the eigenf unctions explicitly |40)| , [41)|1 , we then give a 



proof that the RHO solution indeed provides the exact spectrum of the models. In the final 
part of this section we briefly mention an application of the model to conductance-oscillation 
phenomena in narrow channels ||66)| . |67)|J40)| . 



6.1. Renormalized-harmonic- os dilator solution 

6.1.1. Calogero-Sutherland model 

To depict the essence of the idea, we begin with the Calogero-Sutherland model which 
is given by interacting spinless fermions (or bosons) confined by harmonic potential ITjiPUl , 



1" d 2 , l» 2 2 , _ A(A + 1) nnA . 

where the interaction parameter is assumed to be A > 0. Here, we recall the key property 
common to 1/r 2 models, i.e. the interaction gives rise to the repulsion among energy levels, 
and enlarges the spacing of quantum numbers uniformly. As seen in the previous sections, 
the repulsion effect can be formulated by introducing the step-wise phase shift of the two- 
body S"-matrix in the periodic case. This is the heart of the ABA in which the renormalized 
quantities correspond to the rapidities in the ABA equations. 

For the case of harmonic confinement, the repulsion of energy levels can be taken into 



account by the renormalization of the quantum numbers for oscillators | j38)|| . In the RHO 
approach, all the interaction effects are conjectured to be incorporated into the renormalized 
quantum number, rij, and the energy is given in the expression for free harmonic oscillators, 
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N 1 

E = uY J {n + -). (105) 

3=1 

We note that the renormalized quantum numbers should be related to the conserved charges 
which ensure the integrability of the model, although the explicit relation has not been 
derived yet. As is the case for the periodic case, the repulsion of energy levels are described 
by introducing the step function. Hence the quantum number rij is to be determined by the 
equation ||38)|| , 

N 

rij =Ij + X^2 9{nj - ni), (106) 
i=i 

where Ij(= 0, 1, 2, ■ ■ •) is the bare quantum number, and the step function is introduced such 
that 6(x) = 1 for x > and 8(x) =0 for x < 0. Consequently we obtain the energy for the 
Calogero-Sutherland model as, 

1 N 1 

E = u,[-\N(N-l)+'£ l {I j + -)]. (107) 

Z 3=1 Z 

The ground state is given by the successive quantum numbers Ij = 0, 1, ■ • • , N — 1, and 
particle-hole excitations are described by changing the quantum numbers Ij from those for 
the ground state. Remarkably enough, the above results deduced from the RHO reproduce 



1 


2) 


29 


,34 



6.1.2. SU(u) confined model 



The RHO solution can be applied to more general multicomponent models ||38)|| . As 
an example, let us study the SU(i/) electron model with harmonic confinement, which is a 
variant of the SU(z/) Sutherland model with periodic boundary conditions discussed in the 
previous sections. The Hamiltonian reads | |34)| , ^5)| , |4irj| . 

^ i=i ux i * i=i j>i y^j x i) 

with the spin-exchange operator P^ 13 of two particles (a, (3 = 1,2, •••,//). As usual, to 
complete the diagonalization for the SU(z/) model, it is necessary to introduce the set of 
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renormalized quantum numbers (a = 1, 2, • • • , u) which satisfy the nested algebraic 



equations [[38 )|| 



M 2 Ah 



f = ~ E K*? ~ n&) + A E ~ (109) 



M a M Q _! M a +1 

E - n ( r ] ) + = E « - ^) + E « - (no) 

« 3 

for 2 < a < v. Here the bare quantum numbers 1^ are the non-negative integers (= 0, 1, ■ ■ •) 
which specify v kinds of elementary excitations. In the above equations, the quantity M a = 
Y?8=a was introduced, where Np is the number of electrons with j3 spin (Mi = N = 

The energy is written in the expression for harmonic oscillators, E = uj J^f^i 11 ^ + 1/2), 
in the RHO method. By substituting the nested equations to this formula, and then iterating 
the substitutions, we obtain the final expression for the energy as 138)1 , 



E = uj[-XN(N - 1) + - -M a (M a - 1) + E li a) ))- (HI) 



1x 2 



3=1 



The ground state is described by the successive non-negative quantum numbers 1^ = 
0, 1, ■ • ■ , M a — 1, resulting in the ground-state energy for SU(i>) singlet, 

E^^-XNiN-^ + ^N^}. (112) 

Z Z Q = l 

We note that all the interaction effects are incorporated via the first term of (112), and 
any effects of the interaction do not show up so far as the number of electrons are fixed. 
Hence, the excitation spectrum is described by free oscillators in case of the fixed number of 
electrons. We should note that this does not mean the system to be out of free oscillators, 
as is clearly seen from the level-repulsion effects in the RHO equations. 



6.1.3. hierarchical models 

It is straightforward to apply the RHO solution to the hierarchical family of the confined 
models analogous to (76) [|42)|] , 
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1 A 9 2 1 ^, 22 ^ V a p(V a p + P°f) 

+ L (t-t)i ' (113) 

^ j=i ^ i=i j>j l x j X «J 

where the interaction parameters are the same as those of (77). Applying the RHO tech- 
niques describe above, we obtain the energy as, 

E/ = ujJ2 + ~(\ a - l)M a (M a - 1) + E ifl (H4) 

a=l / / j=l 

The ground-state energy then takes the form, 

E o = % Et^a + A Q M a (M a - 1)]. (115) 

z a=l 

We can classify the excitation spectrum in the matrix formula, which leads to the expression 
similar to (92) in which we should replace (2nv/L) by uo. The detail of this part will be 
reported elsewhere ||42) |. 



6.2. Construction of eigenfunctions 

We have studied so far the energy spectrum making use of the RHO solution which is 
deduced from the level-repulsion effects of 1/r 2 interaction. Here we give the microscopic 
foundation of the RHO solution by constructing the eigenfunctions explicitly. We give a 
brief review of the results of Vacek et al. to construct the eigenfunctions of the ground-state 
i0] \ as well as the excited states [ 41] 



6.2.1. ground state 

We start by writing down a general form of Jastrow wavefunction which is to be the 
exact eigenstate of the family of 1/r 2 models. It consists of two parts, 

&(xiai, . . . , x N a N ) = Y[ \xj — Xi| A ^ (^io;i, . . . , x N a N ), (H6) 

j>i 

where the first one is the Jastrow factor, and the second, is the eigenfunction for the 
noninteracting SU(v) electron model (a = 1, 2, • • • , v). For the above wavefunction, there- 
fore, all the interaction effects are assumed to be taken into account solely by the Jastrow 
factor \xj — Xi\ x . We should like to emphasize that this expression is quite general for the 
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1/r 2 models. For example, the wavefunction for the SU(z/) Sutherland model takes this form 
(99). We will see that this is also the case for excited states. 

For the confined model of (108), the noninteracting wavefunctions is given by that for 
free electrons in the harmonic potential. Hence, we can expect that the ground-state wave- 
function should take the form |34)| , f40)|| , 



n 



i-$gn(oij - a,-, 



N 



10 



Il cx P(-o 



X, 



117) 



i=l 



It has been shown that this wavefunction is indeed the exact eigenfunction of the confined 
model (108) |40 )| . Applying the kinetic term and the potential term on the wavefunction, 
we obtain 



1 N 

~E 



d 2 



Ox; 



+ uj 2 x 2 



00 



\N(N-l) + J2N 2 



u. 



118) 



where the term u is given by 



u 



^ A(A - 1) ^ ^ 2\5 akae 



k<£ 



(X k - Xi) 2 



(x k - Xt)'' 



+ E 



i^k^e 



(Xi - X k )(Xi - x e )' 



119) 



On the other hand, the action of the interaction term yields HO 



£ A(A + Pff) 



j>i 



( Xj -Xi) 2 



^ A(A - 1) ^ ^ 2X5 akae 



k<£ 



[Xk - 



+ E 



A 



k<e [X k - Xi) 



i-n 

ij=U 



Xi — Xi \ "° 
Xi X k J 



[1-S akae ). (120) 



A remarkable point is that multiparticle terms in the above expressions (119) and (120) 
coming from the kinetic energy and the interaction energy cancel each other completely 



10) | . Consequently, it has been proven that the Jastrow wavefunction (117) is the exact 



eigenfunction of the Hamiltonian (108), and the corresponding energy is given by the ex- 
pression (115). The detail of the calculation is given in ref. [ID]]. One can see evidence that 
the above eigenstate indeed corresponds to the ground state, though it is not easy to give a 
rigorous proof for it. For instance, in the case of spinless fermions, the wavefunction (117) 
reduces to the exact ground-state wavefunction [p^], |2]] , p7)[] . Also, in the limit of uj — > 0, 
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we can show that (117) is the ground-state wavefunction. Based on these observations we 
believe that the eigenfunction generally describes the exact ground state of the Hamiltonian. 

6.2.2. excited states 

Concerning the eigenf unctions for excited states, we can also use a general form of the 
Jastrow wavefunction (116). In this case, it is necessary to introduce the exited wavefunc- 
tions of free electrons for \Po. Therefore the following Jastrow wavefunction should be a 
candidate for the exact excited state ||41)|| , 



j>i 



I A 



N 1 , 

1 „ , rmo 



<><y+...-\-H> V- / ' = 1 ^ 



*g, (121) 



where H m is the Hermite polynomials which can produce the excited states systematically. 
It has been shown explicitly in [|41 )|] that this wavefunction gives the exact excited state 
of the Hamiltonian (108) with the corresponding energy obtained by the RHO method 
(114). Therefore, the microscopic derivation of the eigenfunctions establishes that the RHO 
approach proposed here gives the exact solution to the confined 1/r 2 systems. We conclude 
this subsection by stressing that the general form of the wavefunction (116) may be quite 
helpful for constructing the excited states generally for the family of the 1/r 2 systems. 



6.3. Application to conductance oscillations 

Here, we briefly mention some attempts to apply the confined 1/r 2 model to the transport 
of electrons through a narrow channel of the semiconductor nanostructure. We first note 
that it is crucial to take into account the effects of mutual electron interactions to explain 



experiments of conductance oscillations ||68)| , |69)|| . Also, it is known experimentally that one- 



dimensional electrons are confined to a finite segment by impurities or constrictions [ |69)| 



Hence it may be necessary to introduce a interacting electron model with certain confining 
potential. Motivated by the above experiments, the 1/r 2 models with harmonic confinement 
have been applied to conductance oscillations |66)||37)|| , and several characteristic properties 
have been explained. Subsequently, the effect of the internal spin degrees of freedom was 
taken into account correctly [p9)|j40)| , and it was demonstrated that there can be two kinds 
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of periods in the conductance oscillations. We briefly summarize the results of refs. [|39)| , |10) . 



We will be concerned with the SU(2) case of (108) which corresponds to the ordinary electron 
systems. 

In order to consider the transport in narrow channels, a weak coupling is introduced be- 
tween the segment given by the model Hamiltonian (108) and two reservoirs. Let us suppose 
that the conductance is controlled by a resonant tunneling between the one- dimensional 
segment and the reservoirs. Therefore a peak in conductance oscillations occurs when 
the chemical potential of the reservoirs satisfy the relation, fi(N) = E g (N + 1) — E g (N). 
The spacing 5 of two successive peaks in the conductance oscillations is then given by 
5(N) = fi(N + 1) — fi(N). Using the formula (115) for the ground state energy for the SU(2) 
case (electrons), one can see that there appear two independent periods of the conductance 
oscillations j39]jp| , 



5 1 = lu\, 5 2 = uj(\+1), (122) 

reflecting the exchange effect due to the internal spin degrees of freedom. We stress that 
this result improves those previously obtained by Tewari ||66)|| and also by Johnson and 



Payne ||67)|| , who employed the spinless fermion models and concluded a single period for 



the conductance oscillations for any strength of the interaction. 

For the parameters employed in ref. ||67)|| , two periods become 5\ = 7.5u and 5 2 = 
8.5a; in the present model. Hence, the correction due to the exchange effect is small for 
these parameters (strong correlation regime). The exchange effect, however, becomes more 
conspicuous when the interaction becomes weaker (smaller A), and then two periods becomes 
more distinct from each other. In the weak-coupling limit (A — > 0), the present model 
reproduces the results for free electrons in the harmonic well: 5\ — > 0, 5 2 — > ou. It may be 
interesting to study experimentally whether two periods due to the exchange effect can be 
observed in the conductance oscillations. 



7. Summary 

40 



We have reviewed our recent works on the quantum 1/r 2 models. We have successfully 
applied the ABA solution to the multicomponent models such as the SU(i') spin chain, the 
Osp(z/, 1) supersymmetric t- J model, and the hierarchical models related to the FQHE. Ap- 
plying CFT techniques to the ABA equations, we have studied the critical behavior of this 
class of integral models. As for the confined models with harmonic potential, the RHO solu- 
tion has been proposed, which enables us to construct the energy spectrum systematically. 
We have then proven, by explicitly constructing the eigenfunctions, that the RHO indeed 
gives the exact solution to the family of the confined 1/r 2 models. The results have been 
used to discuss the exchange-correlation effects on the conductance oscillations in narrow 
channels. 
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